Flavoring Astrophysical Neutrinos: Flavor Ratios Depend on Energy 
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Electromagnetic (and adiabatic) energy losses of tt's and /i's modify the flavor ratio (measured at 
Earth) of neutrinos produced by n decay in astrophysical sources, : : , from 1 : 1 : 1 at 
low energy to 1 : 1.8 : 1.8 at high energy. The transition occurs over 1-2 decades of v energy, and 
is correlated with a modification of the neutrino spectrum. For 7-ray bursts, e.g., the transition 
is expected at ~ 100 TeV, and may be detected by km-scale v telescopes. Measurements of the 
transition energy and energy- width will provide unique probes of the physics of the sources, vr 
and jj, energy losses also affect the ratio of flux to total 1/ flux, which may be measured at the 
W-resonance (6.3 PeV): It is modified from 1/6 (1/15) at low energy to 1/9 (practically 0) at high 
energy for neutrinos produced in pp {jrf) interactions. 

PACS numbers: 14.60.Pq, 96.40.Tv, 98.70.Rz, 98.70.Sa 



The existence of extra Galactic high energy neu- 
trino sources (for reviews see Q ) is implied by ob- 
servations of ultra-high energy cosmic rays with 
energies of e > 10^^ eV. The acceleration of par- 
ticles in extra-Galactic sources is expected to pro- 
duce a cosmic-ray spectrum extending over many 
decades of energy (although the observed spec- 
trums is likely dominated by Galactic sources be- 
low ~ 10^^ eV), leading to the production of a 
wide energy spectrum of extra-Galactic high en- 
ergy neutrinos. Possible sources of high energy 
neutrinos include 7-ray bursts (GRBs) and active 
galactic nuclei (AGNs). High energy (> 0.1 TeV) 
neutrino telescopes are being constructed in order 
to detect extra Galactic neutrinos (for review see 
0). Their detection will allow one to identify the 
high energy cosmic ray sources, and to probe their 
physics. It may also provide information on fun- 
damental neutrino properties. 

High energy neutrinos are believed to be pro- 
duced in astrophysical sources mainly through the 
decay of charged pions, 7r+ — > /i"'" -I- i/^ — > e"*" -I- 

+ Ve + v^j_ or TT" /i" + Vf^^ er + Vf, + + 
produced in interactions of high energy protons 
with "target" photons (^7) or nuclcons {pp, pn)^. 
The ratio of the fluxes of neutrinos of different fla- 
vors is therefore expected to be, at the source, 
: $° : $0^ = 1 : 2 : ($^, stands for the 
combined flux of vi and Di). Neutrino oscillations 
then lead to an observed flux ratio on Earth of 
<i>y^ : : = 1 : 1 : 1 0. For neutrinos pro- 
duced in pp (pn) collisions, where both 7r'*"'s and 
7r~'s are produced, the ratio of 9^ flux to total v 



v-r^s may also bo produced, by production of charmed 
mesons, but the higher energy threshold and lower cross- 
section for charmed meson production typically imply a 
negligible Vt fraction, e.g. |3| 



flux is 1/6, while for neutrinos produced in p^ col- 
lisions, where only 7r'*"'s arc produced, the ratio is 

1/15 Ulii. 

Km-scale optical Cerenkov neutrino telescopes, 
such as IceGube, are expected to be capable of 
discriminating between neutrino flavors [3| . More- 
over, although these detectors can not distinguish 
between neutrinos and anti-neutrinos at most en- 
ergies, they may allow one to determine the ra- 
tio of Ve and fluxes by detecting v^s at the 
VF-resonance {Vf.e~ W~ anything) around 
6.3 PeV (see, e.g., for discussion of resonant 
9^ detection in IceCube). It has been pointed out 
0, Q that measurements of the flavor composition 
of astrophysical high energy neutrinos may enable 
one to probe new physics, by searching for devi- 
ations from the standard flavor ratio 1:1:1. 
Measurements of the Ve to z/g flux ratio may al- 
low one to probe the physics of the sources [5j , by 
discriminating between the two primary modes of 
pion production, p7 and pp collisions. 

Although a 1 : 1 : 1 flavor ratio appears to be 
a robust prediction of models where neutrinos are 
produced by pion decay, we point out in this let- 
ter that energy dependence of the flavor ratio is a 
generic feature of models of high energy astrophys- 
ical neutrino sources. Pions are typically produced 
in environments where they may suffer significant 
energy losses prior to decay, due to interaction 
with radiation and magnetic fields Since the 

pion life time is shorter than the muon's, at suffi- 
ciently high energy the probability for pion decay 
prior to significant energy loss is higher than the 
corresponding probability for muon decay. This 
leads to suppression at high energy of the relative 
contribution of muon decay to the neutrino flux. 
The flavor ratio is modified to : 1 : at the 
source, similar to that of atmospheric neutrinos at 
high energies where the muons do not decay (see 
e.g. [13), implying 1 : 1.8 : 1.8 ratio on Earth 



The energy dependence of the neutrino flavor 
content provides a unique probe of the sources. On 
the other hand, it comphcates attempts to study 
new physics based on measuring deviations from 
^i^c ■ ^I'ti =1:1- Furthermore, if muon energy 
losses become important at < 1 PeV, it would af- 
fect the flux of De near the W resonance, render- 
ing the suggestion to probe neutrino mixin g an gles 
with neutrinos around the W resonance im- 
practical, and making the discrimination between 
P7 and pp collisions more difficult (due to reduc- 
tion of : $*°*'^'). 

We first derive below approximate analytic ex- 
pressions describing the energy dependence of the 
neutrino flux, flavor ratio and anti-particle con- 
tent, for sources that produce pions with a power- 
law energy distribution [a differential number flux 
$,r(e7r) oc e^^], and assuming the charged par- 
ticle energy loss rate to be proportional to a 
power of the particle energy, e oc — e". Such 
energy dependence is expected for synchrotron 
and inverse-Compton emission (below the Klein- 
Nishina regime), in which case e oc — e^, and for 
"adiabatic" energy loss (energy loss due to ex- 
pansion of the plasma in which the particles are 
confined), in which case e oc — e^. We then dis- 
cuss several specific models of high-energy neu- 
trino sources, and the general implications of our 
results. 

The energy dependence of the flavor ratio. We 
consider neutrinos from astrophysical sources pro- 
duced by the decay of charged pions. The decay 
of the pion and the subsequent decay of the muon 
lead to a flavor ratio of <f>° : $° : $° =1:2:0. 

t^fi 

We note that in the production of a charged pion, 
the high energy proton may be converted to a high 
energy neutron (e.g. p-l-7 — > n+7r+). The neutron 
may escape the source and decay, producing an 
additional Pg. However, in this decay the neutrino 
carries only a small fraction, ^ 10^'^, of the origi- 
nal proton energy, comparable to (7ti„ — mp)/mm 
much lower than the energy of neutrinos produced 
by the pion decay. Since in most astrophysical 
sources the number of emitted neutrinos drops 
rapidly with energy, the contribution of neutron 
decay to the neutrino flux would generally be small 
and is therefore neglected here. 

The flavor ratio $° : $° : $0 =1:2:0 
is modified at high energy, where the life time 
of the muons becomes comparable to, or smaller 
than, the time for significant electromagnetic (or 
adiabatic) energy loss. The particle life time, 
T'a:, decay wcrc X stands for TT^ or /i*, is propor- 
tional to tx- Thus, the ratio of cooling time, 
'nE,cooi = to life time is rapidly decreas- 

ing with energy, Tx,coo\lTxAcca,y « e;^"- We denote 



the energy at which Tx,coo\ = Tx^Accay by eo,x- The 
1/e" dependence of the cooling and decay time ra- 
tio implies that eo,7r/eo,^j is approximately given by 
(ro,^/To,;x)-'/" ~ 102/", where tq,. = 2.6 x 10-^ s 



and To , 



2.2 X 10 s are the pion and muon 



rest-frame life times. 

Wc consider below neutrinos produced by the 
decay of 7r"^'s. The results for neutrinos from ■k~ 
decay are simply given by replacing each particle 
with its anti-particle. For pion decay, each of the 
four final leptons carry approximately 1/4 of the 
pion energy. The (differential number) flux of I'^'s 
of energy produced by 7r+ decay is therefore 
approximately given by 

/•OO 

<(e,) = -9, / de,$,(eOPx(e„4e,). (1) 



Here, $7r(ei) is number flux (per unit energy) of 
pions produced by the source with energy e^, and 
P7r(ei,4e^) is the probability that a pion produced 
with energy would decay before its energy drops 
below 4e^. 

For e oc e" with n > 0, P(ej,e/) = 1 - 
exp[-ejf(ej" - e.^")/"-]- Assuming $^(e,) oc e,^'' 
with A: > 1, we find that the neutrino flux is sup- 
pressed due to energy loss by a factor 

e(^.) 



(no loss) 




' e 



k - 1 



(2) 



Here s = (eo,7r/4e,y)"/n, 7(0, z) is the lower incom- 
plete gamma function, and $q'^° is the flux 
that would have been obtained had energy losses 
been negligible. 

Similarly, for z^'s produced by the decay of /i+ 
we have 



de,, 



X F^(e^,3e,)9,^F,(e„4e^/3), (3) 



and 



(no loss) 
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(4) 



-q s^l- 
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-q s 



where q = 4eo,^/3eo,7r- 

The flavor content given in eqs. (|2I4|I is modi- 
fied as the i/^s propagate from the source to Earth. 
For propagation over cosmological distances d, 
d ^ ftce/Am^c^ and the observed fluxes <l>j,^ 
{a ~ e,fi,T) are related to the production fluxes, 
$0 , by (see e.g. 



P i 



(5) 
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(a)v flavors at the source (b)v flavors on Earth 
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FIG. 1: Flavor and anti-particle content of the flux of astrophysical neutrinos produced by pion decay, for 
$,r oc and oc e^. Figures (a-b) present the energy fluxes in different flavors, e^^^v^ (normalized to e^^^i^e 
at low energy), t&^j stands for the combined flux of ui and vi, and these plots are therefore valid for neutrinos 
produced by any combination of tt^ and tt~ decay. Figure (c) presents the ratio between $^^{^^) and (solid 
line), with 90% CL lines of Vfj, (dashed) and (dotted) fluxes. Figure (d) presents the ratio of to total v flux 
on Earth, solid (dashed) line for neutrinos produced by P7 (pp) interactions. 



Here, Uai is the neutrino mixing matrix, li^a) = 
J2i=i^aiWi) where Vi {i = 1,2,3) are the mass 
eigen-states. 

Uai can be written as a function of three mix- 
ing angles 6*12, 6*23, 6'i3and a Dirac phase S. The 
best fit for (5 = 113 is: 6*12 = 34° ± 2.5°, 
6*23 = 45° ± 6°, 6*13 = 0° ± 8° (90% confidence 
level). The observed neutrino flux ratio is : 
: = 1 : 1.00 ± 0.15 : 1.00 ± 0.15 for ini- 
tial flux ratio : : 1 : 2 : 0, and 
: : = 1 : "1.8^0.70 : 1.8 ± 0.47 for 
: : = : 1 : (90% confidence 
level). Here, we have obtained an approximate 
estimate of the error bars by assuming the errors 
in different mixing angles are not correlated, i.e. 
A($.J<i>.J = {E,,^0^,[deJ'P.J^.J?V/'. 

The results of applying eq. jSJ to the initial flux 
ratios given by eqs. (0 0)) are shown in fig. Qfa- 
c), for the case of 1/e^ and oc e^. As 

expected, the muon (pion) decay neutrino flux is 
suppressed above ~ eo./j/3 (e,y ^ eo.7r/4) by a 
factor oc e"" = e~^, reflecting the fact that at high 
energy the probability for decay prior to significant 
energy loss is oc e^" = e^^, and the flavor ratio 



transition takes place over an energy range of ~ 
eo.^/eo.A. - (ro,./To,^)-i/" ~ 10^/" ^ 10. 

The energy at which a flavor ratio transition 
takes place (from 1 : 1 : 1 to 1 : 1.8 : 1.8), and 
the width of the transition, provide unique han- 
dles on the properties of the source. £0,^1 ni^'Y 
constrain, e.g., the (radiation and magnetic field) 
energy density at the source, and the width of the 
transition may discriminate between adiabatic and 
electromagnetic energy losses. It should be kept in 
mind, however, that while the flavor ratio transi- 
tion is expected to be a generic feature of high en- 
ergy neutrinos produced in astrophysical sources, 
the behavior in the transition region may be more 
complicated that described by eqs. which 
were obtained under idealized assumptions. For 
example, neutrinos may be emitted from different 
regions (in a single source) with different values 
of eo^fj_ (thus "smearing" the transition), and sup- 
pression of the inverse- Compton scattering cross 
section in the Klein-Nishina regime may lead to 
deviations from a simple power-law energy depen- 
dence of the energy loss rate. 

Neutrinos from specific sources. GRBs are pos- 
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sible sources of high energy neutrinos. The 7- 
rays are beheved to be produced by the dissipa- 
tion of the kinetic energy of a highly-relativistic 
wind. Neutrinos that are expected to be produced 
in the same region where the 7-rays are produced 
have characteristic energies > 100 TeV Q- The 
pions and muons cool by synchrotron radiation, 
K , and the muon coohng energy in the 7-ray 
production region is 0, 0| 

eo., « lO^li^TeV. (6) 

Here T = 10^'^r2.5 is the wind Lorentz factor, L = 
lO^^Lsa erg/s is the kinetic energy luminosity of 
the wind (assuming spherically symmetric wind), 
and At = 10~^At_3 s is the observed variability 
time scale of the 7-ray signal. 

If GRBs are associated, as commonly believed, 
with collapses of massive stars, neutrinos of lower 
energy, > 1 TeV, may be emitted in the early 
stage of GRB evolution, when the relativistic wind 
penetrates the progenitor star [l^ . The high en- 
ergy density of radiation implies, in this scenario, 
£0,^ < 1 TeV due to inverse-Compton losses, and 
the Klcin-Nishina suppression of inverse-Compton 
losses at high energy implies eo,7r ^ 1 TeV. In this 
case, an observed flavor ratio of 1 : 1.8 : 1.8 would 
therefore be expected at all energies (> 1 TeV). 

Measurements of the energy dependence of the 
neutrino flavor ratio would therefore provide con- 
straints on the physical parameters of the source, 
and may allow one to discriminate between differ- 
ent scenarios for neutrino production in GRBs. 

Eq. (jnj holds also for neutrinos produced in 
blazar AGN jets For the parameters charac- 
terizing these objects, F ~ 10, L ~ lO''^ erg/s and 
At ~ lO'' s, we have eo,,, 4 x lO*' TeV, imply- 
ing a flavor transition at ^ 10^ TeV. When the 
dominant cooling process is adiabatic cooling, a 
similar cooling energy is obtained. However, at 
these energies, the number of neutrinos detected 
from AGNs may be too small to allow detection of 
the transition. 

Implications for the flux of T>e ■ The fraction of Pg 
flux out of the total flux decreases at high energy 
due to muon cooling. For neutrinos produced in 
■pp interactions, where tt+'s and tt^'s are produced 
at roughly equal numbers, a flavor ratio of 1 : 1 : 
1 is obtained at low energy for both and i/q., 
implying <I>p^ : (Jto**^' = 1:6. At high energy, 
the flavor ratio changes to 1 : 1.8 : 1.8, implying 
^j/e '■ = 1:9. For neutrinos produced in 

P7 interactions, where only tt+'s are created, <l?p^ : 
(j, total _ 2 : 15 at low energy. At high energy, 
the muon energy losses suppress the production of 



anti-neutrinos, resulting in a strong suppression of 
'^p^/^lf^''^ This is illustrated in fig. Eld), where 
the energy dependence of the v^, to Vf^ flux ratio is 
plotted for <I>7r oc 1/eJ and kx e^. 

If the transition is below 1 PeV, as expected for 
GRBs, it will reduce the flux of i^e near the W res- 
onance, making the detection of i^e, and hence the 
discrimination between p7 and pp collisions, more 
difficult. It should be pointed out in this context, 
that significant production of tt^'s may occur in 
sources where pions are produced only through in- 
teractions of nucleons with photons, if the photo- 
production optical depth is large enough, so that 
the neutron produced in a p7 — > n + 7r+ inter- 
action is likely to interact with a photon before 
escaping the source (71.7 p + tt^). Thus, a dis- 
crimination between photo-production and inelas- 
tic nuclear collision sources based on the observed 
^pc/^l"*'^' ratio may not be straight forward. 
Conclusions. We have shown that a flavor ratio 
transition, from : : (f>^^ = 1 : 1 : 1 at low 
energy to 1 : 1.8 : 1.8 at high energy, is expected 
to be a generic feature of high energy neutrinos 
produced in astrophysical sources. The location 
and energy width of the transition provide unique 
handles on the properties of the source, and may 
allow one to discriminate between different scenar- 
ios for neutrino production. The modified flavor 
ratio affects the experimental upper limits on the 
total neutrino flux, which are commonly obtained 
assuming a 1 : 1 : 1 ratio (e.g. 0). It also changes 
the ratio of De to total ly flux from 1/6 (1/15) at 
low energy to 1/9 (practically 0) at high energy 
for neutrinos produced in pp {p"f) interactions. 

We have neglected in the current analysis the 
possibility of matter oscillations, which are typi- 
cally unimportant for the high energy sources un- 
der consideration (e.g. 01). However, it should be 
kept in mind that matter oscillations may, in gen- 
eral, play a role in modifying the predicted flavor 
ratio. For example, the flavor ratio of neutrinos 
produced by a GRB jet buried in a blue super 
giant would be modified by Mikheyev-Smirnov- 
Wolfenstein resonance to 1 : 1.5 : 1.5 outside the 
source and 1.2 : 1.4 : 1.4 on Earth. 
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